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We apply the light-front reduction of the Bethe-Salpeter equation to matrix elements of the 
electromagnetic current between bound states. Using a simple (l + l)-dimensional model to calculate 
form factors, we focus on two cases. In one case, the interaction is dominated by a term instantaneous 
in light-cone time. Here effects of higher Fock states are negligible and the form factor can be 
effectively expressed using non-wave function vertices and crossed interactions. If the interaction 
is not instantaneous, non-wave function vertices are replaced by contributions from higher Fock 
states. These higher Fock components arise from the covariant formalism via the energy poles of 
the Bethe-Salpeter vertex and the electromagnetic vertex. The replacement of non-wave function 
vertices in time-ordered perturbation theory is a theorem which directly extends to generalized 
parton distributions, e.g., in (3 + 1) dimensions. 
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I. INTRODUCTION 

More than a half century ago, Dirac's paper on the 
forms of relativistic dynamics Q] introduced the front- 
form Hamiltonian approach. Applications to quantum 
mechanics and field theory were overlooked at the time 
due to the appearance of covariant perturbation theory. 
The reemergence of front-form dynamics was largely mo- 
tivated by simplicity as well as physicality. The light- 
front approach has the largest stability group || of any 
Hamiltonian theory. Today the physical connection to 
light- front dynamics is transparent: hard scattering pro- 
cesses probe a light-cone correlation of the fields. Not 
surprisingly, then, many perturbative QCD applications 
can be treated on the light front, see e.g. ||. Outside this 
realm, physics on the light cone has been extensively de- 
veloped for non-perturbative QCD Q] as well as applied 
to nuclear physics J|, |J. 

This paper concerns current matrix elements between 
bound states of two particles in the light-front formal- 
ism. We approach the topic, however, from covariant 
perturbation theory. As demonstrated by the tremen- 
dous undertaking of [Q, one can derive light-front per- 
turbation theory for scattering states by projecting co- 
variant perturbation theory onto the light cone, thereby 
demonstrating their equivalence including the delicate 
issue of renormalization. As to the issue of light-front 
bound states, a reduction scheme for the Bethe-Salpeter 
equation recently appeared Q that produces a kernel cal- 
culated in light-front perturbation theory. For the pur- 
pose of simplicity, we consider only bound states of two 
scalars interacting via the exchange of a massive scalar 
in the (1 + l)-dimcnsional ladder model. This work su- 
persedes our original investigation |5). 

Our main consideration is to extend the reduction to 
current matrix elements to investigate valence and non- 
valence contributions in the light-front reduced, Bcthc- 
Salpeter formalism. Thus we calculate our model's form 
factor; moreover, in (1 + 1) dimensions we cannot choose 



a frame of reference where Z-graphs vanish. This enables 
us to completely investigate their contribution, which in 
(3 + 1) dimensions has a variety of applications such as 
generalized parton distributions JhJ. These applications 
are pursued elsewhere [ pd| . Z-graph contributions haunt 
light-front dynamics since non-valence properties of the 
bound state are involved, so that valence wave function 
models cannot be utilized directly. On one hand, the 
light-cone Fock representation provides expressions for 
the Z-graph contributions in terms of Fock component 
overlaps which are non-diagonal in particle number [ fTjJ . 
While on the other hand, vertices which cannot be re- 
lated to the valence wave function (coined as non-wave 
function vertices in fl3|| ) appear in the light-front Bethe- 
Salpeter formalism. A variety of ways have been pro- 
posed for dealing with these non-wave function vertices 
fl4| , |l5| , p| [T7[ | . We also note that one can avoid the issue 
by attempting to model the covariant vertex |l8|, , or 
(when possible) estimating the contribution from higher 
Fock states |(J. 

Below we show that non-wave function vertices are 
supplanted by contributions from higher Fock states in 
light-front time-ordered perturbation theory (provided 
the interaction has light-cone time dependence). In 
essence contributions from non-wave function vertices are 
reducible and should only be used when the interaction 
is (or is approximately) instantaneous 1 . This constitutes 
a replacement theorem for non-wave function vertices 
which trivially extends to (3+1) dimensions. 

The organization of the paper is as follows. First in sec- 
tion [H] we present the issue of non-wave function vertices 
and energy poles of the Bethe-Salpeter vertex focusing 
on the light-front constituent quark model as an cxam- 



We shall often refer to interactions and vertices merely as in- 
stantaneous if they are independent of light-cone time, or equiv- 
alently light- cone energy. 



pie. Non-wave function vertices are required to express 
the form factor. We find the commonly used assumptions 
in quark models necessitate vertices not only without en- 
ergy poles but without energy dependence. Next in sec- 
tion we review the reduction of the Bethe-Salpeter 
equation presented in |§| focusing on the energy poles of 
the vertex. We derive an interpretation of the reduction 
as a procedure for approximating the poles of the ver- 
tex. Additionally in this section we construct the gauge 
invariant current to be used with the reduced formal- 
ism. In section [TV] the ladder model is presented and we 
compare the calculation of the form factor for the model 
using two different paths to the reduction. The compari- 
son allows us to see when non-wave function vertices can 
be efficiently used. Lastly we summarize our findings in 
section 



II. POLES OF THE BETHE-SALPETER 
VERTEX 

To introduce the reader to non-wave function vertices 
and instantaneous approximations, we focus on light- 
front constituent quark models. We start by writing 
down the covariant equation for the meson vertex func- 
tion T. It satisfies a simple Bethe-Salpeter equation H] 
(see Figure [I]): 



/j2 
^JLv(k, P )y BS (p,R), 



(1) 



in which we have defined the Bethe-Salpeter wave func- 
tion f bs as 



V BS (k,R) = G(k,R)r(k,R), 



(2) 



with the two-particle disconnected propagator G(k, R) = 
d(k)d(R — k). For scalars of mass to, the renormalized, 
single-particle propagator d has a Klein-Gordon form 



d(k) 



{k 2 - TO 2 ) [I + (fc 2 - TO 2 )/(fc 2 )] + it 



(3) 



where the residue is i at the physical mass pole and 
the function f(k 2 ) characterizes the renormalized, one- 
particle irreducible self-interactions. To simplify the 



comparison carried out in section [V, we shall ignore 
f(k). Above V is the irreducible two-to-two scattering 
kernel which we shall refer to inelegantly as the interac- 
tion potential. 

Now we imagine the initial conditions of our system 
are specified on the hypersurface x + = 0. We de- 
fine the plus and minus components of any vector by 
x ± = ± x 3 ). Correlations between field opera- 

tors evaluated at equal light-cone time turn up in hard 
processes for which this choice of initial surface is natu- 
ral. In order to project out the initial conditions (wave 
functions, etc.) of our system, we must perform the inte- 
gration over the Fourier conjugate to x + , namely k~ . For 



V 



FIG. 1: Diagrammatic representation of the Bethe-Salpeter 
equation. The blob represents the vertex function V. 



instance, our concern is with the light- front wave function 
defined as the projection of the covariant Bethe-Salpeter 
wave function onto x + = 0, 



ip(x) = 2R + x(l 



dk~ 

, ) / —V BS (k,R) 7 (4) 



with x = k+ /R+. 

Looking at Eq. (Q), in order to project the wave func- 
tion exactly, we must know the analytic structure of 
the bound-state vertex function. If the vertex function 
r(fc, R) had no poles in k~, then our task would be simple: 
the light-front projection of ^bs would pick up contri- 
butions only from the poles of the propagator G(k,R). 
Next we observe from Eq. ([j]), that the k~ dependence 
of the interaction V{k 1 p) must give rise to the k~ poles 
of the vertex function T(k,R). Hence an instantaneous 
interaction gives rise to an instantaneous vertex and a 
simple light-front projection (see e.g. juj). 

On the other hand, constituent quark models often 
assume a less restrictive simplification of the analytic 
structure of the vertex (see e.g. [^9|) in order to per- 
mit the light-cone projection. We shall show that this 
assumption along with the presumed covariance of the 
quark model often implies instantaneous vertices. For 
any momentum p, let us denote the on-shell energy 
Pon = "i 2 /2p + . The propagator G(k,R) has two poles 



on x 

R~ + (k- R)i 



(5) 



Notice that although we work in (I + 1) dimensions, the 
results are trivial to extend to (3+1) dimensions because 
the imaginary parts of poles have precisely the same de- 
pendence on (only) the plus-momenta. This remark ap- 
plies not only to this section but to the rest of this work. 

In constituent quark models T is assumed to have no 
poles in the upper-half /c^-plane for < x < 1. Since the 
light-front wave function is proportional to 9[x(l — x)], 
in light of Eq. (|b|) we further require any poles of T to lie 
in the upper-half plane for x < and in the lower-half 
plane for x > 1. With these restrictions Eq. (^) dictates 
the form of the constituent quark wave function 



tP(x) = D w (x\M 2 )r(k b ,R)0[x{l - x)}, 

where we have defined the Weinberg propagator as 

1 



D w (x\M 2 ) 



M 2 



(6) 



(7) 



x(l—x) 



3 




FIG. 2: Covariant diagram for the electromagnetic form fac- 
tor. Labels a, b and c denote subscripts used for on-shell en- 
ergy poles. 



and used the abbreviation r(fc{,, R) to denote evaluation 
at the pole k~ — k^ appearing in Eq. (jsj) . 

When we calculate the (elastic) electromagnetic form 
factor for these constituent quark models (see Figure]^), 
we are confronted with more poles. Let A + = — C,P + < 
be the plus-component of the momentum transfer be- 
tween initial and final state mesons. The form factor is 
then 



F(t) oc 



(2x - C)r(fe, P)T* (fc + A, P + A) d 2 k 
[fc 2 - to 2 ][(P - fc) 2 - m 2 ][(k + A) 2 - to 2 ] 

(8) 

defined in 



The fc -poles from the propagators are k a , k b 
(§) and fc c 7 = - A- + (fc + A)" - ^ . 

Given this pole structure, the contributions to F(t) are 
proportional to 6[x(l — x)]. In the region x < C, closing 
the contour in the upper-(lower-) half plane will enclose 
possible poles of the final (initial) vertex. As one can see 
from considering the region x > the form factor can be 
determined solely from Res(fc^). In (3 + 1) dimensions, 
where we are free to choose frames in which £ = 0, this 
is the only contribution to the form factor. But if £ ^ 0, 
the additional poles from vertices in the region x < £ are 
required by Lorentz invariance. 

The authors |l5| advocate no modification to the pole 
structure of Eq. (||) due to the vertices. Closing the con- 
tour in the lower-half plane, they pick up the residue at 
k~ without any contributions from poles of the initial- 
state vertex. Such poles cannot lie in the upper-half 
plane for x < £ since the form of tp would be both frame 
dependent and contrary to that of Eq. (^) (the (3 + 1)- 
dimensional version of which is employed by |l5|| ). Thus 
the authors are actually assuming there are no poles of 
the Bethe-Salpeter vertex [|l6| . 

Returning to the definition of the wave function (^) un- 
der the premise of a vertex devoid of poles, we now find 
r(fcf,,P) = r(fc ,P), which we shall call pole symmetry. 
This symmetry is essential for making contact with the 
Drell-Yan formula since for x > £ both initial- and 
final- state vertices may be expressed in terms of wave 
functions. When x < £, however, the final-state vertex 
becomes F*(fc a + A, P + A) which cannot be expressed 



in terms of ip*(x') (where x' = j^) even with pole sym- 
metry. Such a vertex is referred to as a non-wave func- 
tion vertex. Understanding and dealing with such objects 
from the perspective of time-ordered perturbation theory 
is the primary goal of this paper. 

Because these constituent quark models are presumed 
covariant, Eq. (||) converges. This enables us to relate an 
initial-state non-wave function vertex to the final-state 
non-wave function vertex encountered above via inte- 
grating around a circle at infinity. Equating the sum 
of residues Res(fc~) + Res(fc^) + Res(fc~) = 0, we find 

rn*. + A,n (v _,- ) _rfc|) s 



r*(fc b + A,P') 



T(k b ,P) 



) — k a k c . 

(9) 

This relation holds for the class of models for which the 
pole structure of Eq. (||) is not modified by the vertices. 
The ratio structure of the vertices does not allow for a 
common factor fc~ — fc~ in the two terms in Eq. @ . The 
equality then depends on delicate cancellations between 
initial- and final- state vertices, which in general are un- 
related. The philosophy of constituent quark model phe- 
nomenology is to choose the form of ip and hence the 
form of r. Treating T free, equality can only hold if both 
ratios are one. This yields the restriction 



d 



dA- 



■r(fc c ,p) = 0. 



(10) 



But A~ only enters r(fc c , P) through k~~ = k~ . Hence T 
is independent of fc - . In general, of course, the vertex T 
not only has light-front energy dependence but poles as 
well. As we shall see below, the light-front reduction of 
the Bethe-Salpeter equation is a procedure for approxi- 
mating the poles of the vertex function. Moreover when 
applied to current matrix elements, these poles generate 
higher Fock state contributions. 



III. LIGHT-FRONT REDUCTION 

In this section, we review the reduction scheme set 
up in Q since their notation, while useful, is unfamiliar. 
Additionally we resolve a peccadillo intrinsically related 
to non-wave function vertices. Next we give an intuitive 
picture for the reduction and finally construct the gauge 
invariant current for the calculation of form factors on 
the light-front. 

Above we have removed overall momentum-conserving 
delta functions, e.g. our propagator G(fc, R) is the mo- 
mentum space version of G(R) defined by (R'\G\R) = 
(2tt) 2 S ( - 2 \R' - R)G(R). In terms of these fully two- 
dimensional quantities, the Lippmann-Schwinger equa- 
tion for the two-particle transition matrix T appears as 



T = V + VGT. 



(11) 



A pole in the T- matrix (at some P 2 = M 2 , say) corre- 
sponds to a two-particle bound state. Investigation of 



4 



the pole's residue gives the Bethe-Salpeter equation for 
the bound-state vertex T 



VGT. 



(12) 



The Bethe-Salpeter amplitude * is defined to be GY. 
Following ||, we denote quantities able to be rendered 
in position or momentum space with bras and kets. 
We will employ this notation only for quantities that 
have been stripped of their overall momentum-conserving 
delta functions, for example T(k,R) — (k\Tn), where R 
is used as a label for the bound state for which R 2 = M 2 . 



A. The scheme 

To reduce the Lippmann-Schwinger equation to a light- 
front version, we must introduce an auxiliary Green's 
function G in place of G (as in (24|]). Thus we have 



provided that 



T = W + WGT, 



W = V + V(G - G)W. 



(13) 



(14) 



Taking residues of Eq. ( |l3| ) gives us an alternate way to 
express the bound state vertex function 



T = WGT. 



(15) 



To choose a light front reduction, G must inher- 
ently be related to projection onto the initial surface 
x + = 0. For simplicity, we denote the integration 

J^(k-\0(R) = O(R). With this notation, we will 

always work in (1 + l)-dimcnsional momentum space for 



which the only sensible matrix elements of 



O(R) are of 



the form (k+ | 0(R)\ p",p + ). The operator 0{R) is de- 
fined similarly. For a useful reduction scheme (one that 



preserves unitarity), we must have G(R) 



G(R) 



The simplest choice of G that results in time-ordered per- 
turbation theory requires 



G(i?) = G(R) g-\R) G{R), 



(16) 



where the reduced disconnected propagator g(R) is de- 
fined by the matrix elements 



(xR+\ G(R) \yR + ) = (xR + \g(R)\yR+), 



explicitly this forces 



(17) 



(xR + \g{R)\yR + ) = 2ttS(xR + - yR + ) 

2m 



x 0[x(l - x)] 



2R+x(l - x) 



D w (x|i? 2 ). (18) 



The inverse propagator g _1 (R) can then be constructed, 
bearing in mind g^ 1 (R) only exits in the subspace where 
g{R) is non-zero. Explicitly 

{xR + \g- 1 (R)\yR + ) = 2nS{xR + - yR + ) 

x 6[x(l - x)] 2R+X } 1 . - X) Dw' 1 ^ 2 )- (19) 



2m 



This forces 



(xR + \g- 1 {R)g(R)\yR + ) = 2n5{xR + - yR+)9[x(l - x)], 

(20) 

which is unity restricted to the subspace where the oper- 
ators g{R) and g~ 1 (R) are defined. We are more careful 
about this point than the authors [|| since the conse- 
quences of Eq. ( |l9| ) are essential for dealing with instan- 
taneous interactions. Notice G(R) defined in Eq. ( |l6| ) is 
non-zero only for plus-momentum fractions between zero 
and one. 

The reduced transition matrix t(R) is 

t(R) = g-\R) G(R)T(R) G(R) g~\R) (21) 



Taking the residue of Eq. (|2l]) at R 2 = M 2 , gives a ho- 
mogeneous equation for the reduced vertex function 7 



\ lR )=w(R)g(R)\ lR ), (22) 
where the reduced auxiliary kernel is 

w(R) = g-^R) G(R)W(R)G(R) g _1 (i?). (23) 

Given this structure, the reduced kernel w(x,y\R 2 ) = 
(xR+\w(R)\yR+) will always be oc 6[x(l - x)]9[y(l - y)]. 
Moreover the reduced vertex 7( 2; \M 2 ) = (xR + \jr) oc 
9[x(l — x)] as a result of Eq. (p2|). 

From (E2) we can define the light-front wave function 
\iJjr) = g(i?)|7fl), notice this too restricts the momen- 
tum fraction x: ip(x) oc 0[x(l — x)]. By iterating the 
Lippmann-Schwinger equation for T twice, it is possible 
to relate T to t and thereby construct T given t, which 
is clearly not possible from the definition (^l[). Tak- 
ing the residue of this relation between T and t yields 
the reduced-to-covariant conversion between bound-state 
vertex functions, namely 



\T R ) = W(R)G(R) |7«> 



(24) 



Finally, we can manipulate the covariant Bethe-Salpeter 
amplitude into the form 

|*a)= (l+(G(R)-G(R))w(R) ) jG(R)\ | 7fl >, (25) 

which justifies the interpretation of \^Pr) as the light-front 
wave function since \^r) = IV'fl)- 



While all light-front reduction schemes when summed 
to all orders yield the x + = projection of the Bcthc- 
Salpeter equation, the choice of G(i?) in Eq. Jl6| ) gener- 
ates a kernel calculated in light-front time-ordered per- 
turbation theory. Lastly the normalization of the covari- 
ant and reduced wave function is discussed in flllfl. 



B. An interpretation for the reduction 

The heart of our intuition about the light cone lies in 
integrating out the minus-momentum dependence of the 
covariant wave function. So we merely cast the formal re- 
duction in a way which highlights the contributions from 
poles of the vertex function. 

Utilizing Eqs. (15ft and (pj), we can show 



G(R)\T R ) = G(R) | 7fl ). 



(26) 



Thus the appearance of G(R)\Tr) has the form of an 
instantaneous approximation since Eq. ( |2^ ) shows that 
it has no minus-momentum poles besides those of the 
propagator G(R). 

This instantaneous approximation appears in deter- 
mining the light-front wave function. Using Eqs. ( |l4| ) 
and (p5|), we have 

\tPr) = \g(R)V(R) 

OO 

x [(g(R) - G(R))V(R)] J G(R)\ | 7H >. (27) 

3=0 



From truncating the series in G(R) — G(R) at some j = 
n — 1 and using a consistent approximation to Eq. (24), 
we are led to the approximate solution 



\ip( n )) = (GV) n G |r), 



(28) 



after having used Eq. (|26|). Thus at any order n in the 
formal reduction scheme, we have iterated the covariant 
Bethe-Salpeter equation n-times and subsequently made 
an instantaneous approximation via Eq. (|26|). Retain- 
ing the minus-momentum dependence in V to n th order 
allows for an n th -order approximation to the vertex func- 
tion's poles. 



C. Current in the reduced formalism 

Here we extend the formalism presented so far to in- 
clude current matrix elements between bound states. We 
do so in a gauge invariant fashion following pjj] . Our no- 
tation, however, is more in line with the elegant method 
of gauging equations presented in [ p6| . This latter gen- 
eral method extends to bound systems of more than two 
particles. 





FIG. 3: Graphical depiction of the irreducible five-point func- 
tion r M as sum of impulse terms and a gauged interaction. By 
construction, T M is gauge invariant. 



by 



Consider first the full four-point function G^ defined 



G (4) = G + GTG. 



(29) 



For later use, it is important to note that the residue of 
G ( ^(R) at the bound state pole R 2 = M 2 is -i\$!n)(®R\- 
Using the Lippmann-Schwinger equation for T, we can 
show the four-point function satisfies 



G (4) = G + GVG {A) . 



(30) 



To discuss electromagnetic current matrix elements, we 
will need the three-point function d^ where the label i 
denotes particle number. We define an irreducible three- 
point function in the obvious way 



(31) 



Now we need to relate the one-particle electromagnetic 
vertex function to the T matrix. Let denote the elec- 
tromagnetic coupling to the constituent particles (since 

our particles are scalars j M = Since the electro- 

magnetic three-point function is irreducible, we have 



rr = G- 1 G {i) f 



(32) 



and by using the definition of G^ Eq. |29|, we have the 
desired relation 



^ — AH 



f + TGf 



(33) 



Notice the right hand side lacks the particle label i. In 
the first term, the bare coupling acts on the ith particle 
while in the second term the bare coupling does not act 
on the ith particle. For this reason we have dropped the 
label which will always be clear from context. 

In considering two propagating particles' interaction 
with a photon, the above definitions lead us to the im- 
pulse approximation to the current 



i — 1 1 u 2 



d^r$. 



(34) 



Additionally the photon could couple to interacting par- 
ticles. Define a gauged interaction U M topologically by 
attaching a photon to the kernel in all possible places. 
This leads us to the irreducible electromagnetic vertex 
defined as (see Figure ||) 



(35) 





+ 




FIG. 4: Graphical depiction of the five-point function G^ 5 ' M . FIG. 5: Diagrammatic representation of the one-boson ex- 
The irreducible five-point function is the gauge invariant r M . change potential w(x,y\M 2 ) appearing in Eq. (fil]). 



which is gauge invariant by construction. 

Lastly to calculate matrix elements of the current be- 
tween bound states it is useful to define a reducible five- 
point function (see Figure ||) 



G (5) M = G (4) r"G (4) . 



(36) 



Having laid out the necessary facts about electromagnetic 
vertex functions and gauge invariant currents, we can 
now specialize to their matrix elements between bound 
states by taking appropriate residues of Eq. (p6|). The 
form factor is then 

-HP'" + P")F(t) = (*j»|r"(-A)|4>}, (37) 

where P'» = P" + and t = A 2 . 

IV. TWO MODELS 
A. Wave functions 

To say anything less than general, we must know the 
minus-momentum dependence of the interaction. We 
therefore adopt a weakly coupled, one-boson exchange 
model for V (the so-called ladder approximation). Sup- 
posing the boson mass is \i and the coupling constant g, 
we have 



V(k,p) 



(p — k) 2 — fi 2 



(38) 



where the energy pole (with respect to p) of the interac- 
tion is 



Pv 



2(p+-k+) 2(p+-k+) 



(39) 



This interaction is non-local in space-time and hence 
does not have an instantaneous piece. The reduced ker- 
nel Eq. (^3| ) is consequently made up of retarded terms 
(i.e. dependent on the eigenvalue M 2 ) where higher order 
in G — G means more particles propagating at a given in- 
stant of light-cone time (see, e.g. ||). The leading-order 
equation for ip from Eq. fl22fl is the Weinberg equation 




1 w(x,y\M 2 )ip(y) 
2(27r)y(l - y) 



dy. (40) 



with the time-ordered one-boson exchange potential (see 
Figure ph calculated to leading order from Eq. (E3h 



w(x,y\M 2 ) 



-a 



x-y 



■9[x(l ~ x)]9[y(l - y)] 



(x-y)D(x,y\M 2 )-{x~y} , (41) 



where 



L> _1 (x,y|M 2 ) = M 2 - 



x — y 1 — x 



(42) 



We can obtain Eq. (^) most simply b y iter ating the 
Bethe-Salpeter equation once (see section [II B ) and then 
projecting onto the light cone. 

In the limit /j 2 ^> m 2 , M 2 , the interaction becomes 
approximately instantaneous, which suggests we separate 
out an instantaneous piece V Q : 



V(k,p) = V + [V(k,p)-V ], 



where 



V = V (x, y) = 



-g 2 9{x-y) 



g 2 0(y - x) 



E 2 {x — y) — /i 2 E 2 (y — x) — p/ 



(43) 



, (44) 



with E as a parameter to be chosen. Of course other 
choices of V are possible. We choose the above form of 
V for two reasons. First is the form of the instantaneous 
approximation wave function 4>{x). When we write the 
potential as Eq. (fl3|), we expand Eq. (|l4| ) to first order 
in g 2 as 

W = V + (V - V„) + (V - V )(g - g)(V- Vo). (45) 

To zeroth order in G — G and V — V , the equation for cj> 
is 



D w (x|M 2 



1 w (x,y)(t)(y) 
o 2(27r)y(l - y) 



dy, (46) 



where by Eq. (Pq), the reduced instantaneous poten- 
tial is merely w (x,y) = 9[x(l - x)]9[y(l - y)]V a (x,y). 
In the instantaneous limit, fi 2 ^> m 2 ,M 2 solutions to 
Eqs. ( pfj| ) and ( |46| ) coincide — both wave functions ap- 
proach ~ Dw(i,4m 2 ). Secondly, we preserve the con- 
tact interaction limit by excluding from V the factor 
9[x(l - x)]9[y(l - y)]. That is, when {^ 2 ,g 2 } -> oo with 
g 2 1 ' jji 2 fixed (or equivalently E — > 0) we have the contact 



7 



interaction 2 V — > g 2 /^ 2 , which rightly knows nothing 
about the momentum fractions x and y. 

We shall now investigate contributions to the form fac- 
tors for each of the models Eqs. ( p8|) and (|i4|). 



B. Form factors 

From Eq. (|37]) we can calculate the form factors for 
each of the models ( |38| ) and (Q) in the reduced formal- 
ism. Working in perturbation theory, we separate out 
contributions up to first order by using Eq. ( |25|) to first 
order in G — G and T M (jjfj) in the first Born approxima- 
tion. The matrix element = {^ P , |r^(-A)|* P ) then 
appears for a model with some kernel V 



J» « ( lP , | |g(P') (l + V(P')(G(P') - G(P')f 
x (^(-A)^ 1 + V r (-A)G(-A)9^(-A)d" 1 ) 
x (l + (G(P)-G(P))F(P))G(P)| \j p ) 

= (Jlo + SJ t + SJ f + SJ %) + °i v % ( 47 ) 
with the leading-order result 

Jlo = (7P'I \G(P')d»(-A)d^G(P)\ \ 7p ). (48) 
The first-order terms are 

6J? = (>yp,\ \G(P')d»(-A)d^ 

x (g(P)-G(P))v(P)G(P)\ | 7P ) 

Jj; = ( 7 P'| |g(P')^(P / )(g(P / ) - G(P') 

xd^(-A)d- 1 G(P) | 7P ) 

*J£ = ( 7P ,| |g(P')(V(-A)G(-A)£"(-A)) 

xd^P)! | 7P ). (49) 

The labels indicate the intu ition behind the reduction 
scheme (seen in section [II B): the term SJj arises from 



one iteration of the covariant Bethe-Salpeter equation 
for the final-state vertex followed by an instantaneous 
approximation (p6|), 8 J? arises in the same way from 
the initial state and 6 comes from one iteration of the 
Lippmann-Schwinger equation (111]). 




FIG. 6: The Z-graph confronting evaluation of the electro- 
magnetic form factor 



Because the leading-order expression (|48|) is indepen- 
dent of the kernel V, the result will have the same form 
for both models. Using the effective resolution of unity, 
the above expression converts into 



-To — 



d 2 p d 2 k 



(7P'|P + ) 



(2tt) 2 (2tt) 2 

x (p\G(P')d + (-A)d^G(P)\k) <k+ | 7P ). (50) 

Bearing in mind the delta function present in G(i?), we 
have the factor 

( P \d + {-A)\k) = -i(2k + +A + )(2w) 2 S 2 (p-k-A), (51) 

Now define x = k + /P + and A + = —(P + as above and 
denote the reduced vertex (k + | 7 p) = 7 (a;|M 2 ). The 
leading-order contribution is then 



'LO 



d 2 k 



j*(x'\M 2 )d(k 



x (2x - C)d(k)d(P - k)j(x\M 2 ), (52) 



where 



x' = jz^- Eq. ( |5^ ) is quite similar to ([; 
x > C evaluation is straightforward and leads to 



For 



F L o(t) = 



6(x-() f dx 2x-( 



1 - C/2 J 2(2tt) x(l - x)x' 

for the non-instantaneous case. For the instantaneous 
case FLo(t), replace tp with <j). 

On the other hand, we know from section [fj the region 
x < Q contains a non-wave function vertex (see Figure |TJ) . 
However, evaluation of Eq. ( p2|) with reduced vertices is 
subtly different than in section ttj. Quite simply, the term 
j*(x'\M 2 ) = by virtue of Eq. (f22|) because x' < 0. Thus 
there is no contribution at leading order for x < C- 

The first-order terms depend explicitly on the interac- 
tion and will hence be considerably different for each of 
the models. We consider each model separately. 



Non-instantaneous case 



2 This exactly soluble (1 + l)-dimensional light-cone model, in 
which ip(x) oc ~D\fi/(x\M 2 ), has been considered earlier in |27|. 



Evaluating contributions at first order for the non- 
instantaneous interaction Eq. ( |38| ) is complicated by the 
presence of poles in the interaction (c/ Eq. (|39|)). First 
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FIG. 7: Diagrams which contribute to the form factor to first 
order in G — G for the non- instantaneous case (for x > £). 




FIG. 8: The remaining diagrams (characterized by x < £) for 
the electromagnetic form factor at first order in G — G for the 
non-instantaneous case. 



we evaluate the first Born term SJif in Eq. (^jj). After 
careful evaluation of the two minus-momentum integrals, 
we have the contribution to 8 for x > £ 



8J+ = 



(x — C) dxdy (2x — £) 



167r 2 xx'y(l — y)y' 



x D{y',x'\M' 



y-x 



D(y,x\M 2 )^(y), (54) 



where y' = fzj- This contribution corresponds to dia- 
gram A in Figure 0. Additionally using a;" = x/f, we 
have for x < C 



6>(C-x) rfxcfy (2x - Q/C 
167r 2 y(l -y)y'x"(l - x") 
g 2 0(y - x 



i>*(y') 



x D w (x"|t)- 



y-x 



D(y,x\M 2 )ip(y), (55) 



which corresponds to diagram D in Figure 

The initial-state iteration term 6 J? in Eq. ( f49|) is com- 
plicated by the subtraction of the two-particle reducible 
contribution 

-{ lP ,\ \G(P')d^-A)d^ G(P)V(P)G(P)\ | 7P ) 

= -( 7P ,| ^(PO^n-A)^^^)! | 7P ). (56) 

Thus this term merely removes contributions which can 
be reduced into the initial-state wave function. Evalua- 
tion of the two minus-momentum integrals yields a con- 
tribution for x > C to S J 4 + 



8J+ 



(x — C) dxdy (2x — C) 



167r 2 xx'(l — x')y(l — y) 



xD(y>,x'\M 2 / 0{y - x) 
y-x 



D(y,x\M 2 )i;(y), (57) 



which corresponds to diagram B in Figure 0. On the 
other hand, for x < C, we have the non-wave function 
vertex j*(x'\M 2 ) for the final state, which vanishes. Sim- 
ilarly the G(P) term vanishes. Thus there is no contri- 
bution to <5Jf for x < £• 

Finally there is the final-state iteration term 8Jf in 
Eq. ([49|). There are only two types of contributions. For 
x > C, that which can be reduced in to the final-state 
wave function is subtracted by G. The remaining term 
is: 



8J+ = 



6{x-Q dxdy'(2x-() 



167T 2 x(l — x)x'y'(l — y 1 ) 



x D(y',x'\M' 



y-x 



D(y,x\M 2 )^(x), (58) 



where y = £ + (1 — Qy' ■ This corresponds to diagram 
C in Figure 0. For x < £, the subtraction term vanishes 
since x' < for which G(P') = 0. The remaining term 
in 8 J~j gives a contribution 



6(C-x) dxdy'(2x-()/( 



16tt 2 (1 - x)x"(l - x")y'(l - y 1 ) 
g 2 9(y-x) 



x D w {x"\t)- 



y-x 



D(y,x\M^(x), (59) 



which corresponds to diagram E in Figure To sum- 
marize, the non- valence correction to the form factor in 
the non-instantaneous case is 



SF 



1 



NI 



l-C/2 



8 J\ + 8 J% 



Sj£ + SJ£ + 8 J 



(60) 



and there arc no non-wave function terms. 



2. Instantaneous case 

The case of an instantaneous interaction is quite dif- 
ferent due to the absence of light-front energy poles in 
V . Note we are working with Eq. ([|5|) to zeroth order in 
G — G and V — V . The first term in Eq. ( |49| ) we consider, 
is the Born term SJif. The pole structure leads only to a 
contribution for x < Q: 



8J+ 



6>(C - x)0(y - x) dxdy (2x - C)/C 
167r 2 y(l - y)y'x"(l - x") 
xP(y')D w (x"\t)V (x",y")<p(y), 



(61) 

and is depicted by the first diagram in Figure |. The 
interaction along the way to the photon vertex is crossed 
and represents pair production off the quark line (x" > 
1)- 

The next term 8 J? simplifies considerably due to the 
absence of light-cone time dependence in V a . As above, 
the final state vertex restricts ( < x < 1. But then we are 



confronted with a factor 
x > 0. Thus 8 J? = 0. 



G(P) - G(P) 



= since 
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FIG. 9: Diagrams with crossed interactions necessary to cal- 
culate the form factor in the region x < £ for the instanta- 
neous case. 



The last term we must consider is o~Jf, 



have the factor 



in which we 
This is zero 



(G(P')-G(P'))\k + A). 

for x - ( > 0, else G(P') = by virtue of Eq. (Jl6j) and 
Eq. (|li^). Thus we only have a contribution for x < £ 

which is from G(P')\k + A). The expression for this 

contribution is 



e((-x)9(y-x)dxdy'(2x-()/( 
16tt 2 (1 - x) x"(l - x") y'(l - y') 
x cf>*(j')w (y',x')D w (ic"\t)^(x), 



(62) 



and is depicted on the right in Figure The interaction 
again is crossed {x' < 0) and represents pair production. 
With Eqs. ( |6l"| ) and (|6|), we have both bare-coupling 
pieces of the full Born series for the photon vertex (fur- 
ther terms in the series, which result from higher order 
terms in the expansion of W, add interaction blocks to 
each diagram on the quark-antiquark pair's path to an- 
nihilation). In this way we recover the Green's function 
from summing the Born series [|l4| |16] |. Notice also the 
above form of 6 is what one would obtain from ex- 
tending the definition of 7 as a non-wave function vertex 
|jl4j| . For the case of an instantaneous interaction, the 
light-front Bcthe-Salpeter formalism automatically incor- 
porates crossing. 

To summarize, the non-valence contribution to the in- 
stantaneous model's form factor is 



6F I = (5J+ + SJ+)/(l~C/2), 



(63) 



and involves crossed interactions or equivalently non- 
wave function contributions. 



3. Comparison 

Now we compare the form factors for the cases of in- 
stantaneous and non-instantaneous interactions. In the 
instantaneous limit /1 2 3> m 2 ,M 2 , however, we under- 
stand the behavior of the wave functions. The wave func- 



tions Eq. ( |40[ ) and Eq. (46) both become narrowly peaked 
about x = 1/2 in the large /j? limit, cf the behavior 
of Dw(i,4m 2 ). Since we are investigating non- valence 
contributions to form factors, we choose additionally to 
solve for the instantaneous wave function to first order in 



2 / 2 


4 - M 2 /m 2 


0.100 


6.57 x 10" 1 


0.316 


2.72 x 10 _1 


1.00 


6.34 x 10~ 2 


3.16 


8.79 x 10" 3 


10.0 


9.57 x 10~ 4 


31.6 


4.81 x 10 -5 



TABLE I: Numerical solution of the bound-state equation 
Eq. dip] ) for various values of [i 2 . The coupling constant a = 
0.100. 



V — V (see Eq. |45|) to put the wave functions on equal 
footing: i.e. <fi(x) rj ip(x), because the difference V — V 
is presumed small. This has the efficacious consequence 
of producing identical leading-order terms cf Eq. ( p33| ) 
and eliminates the issue of normalization. The optimal 
choice of the instantaneous interaction is not under in- 
vestigation here. So we shall simplify matters further by 
choosing E = in Eq. (|||). 

For a few values of /x 2 , we solve for the wave function 
using Eq. (^0|) . In Table |, we list the values of /i 2 used as 
well as the corresponding eigenvalue M 2 (all for the cou- 
pling constant a = 0.100, where a = g 2 /Amrn A ). We then 
calculate the form factors in the instantaneous and non- 
instantaneous cases. We arbitrarily choose £ = 0.707. In 
(1 + 1) dimensions, this fixes A 2 = -( 2 M 2 /{1 - Q. The 
ratio of these form factors — defined using Eqs. ( p3| , |60| , p3[ ) , 
namely 



(FLo + SF^/iFLo + SFNj) 



(64) 



— is plotted versus log 10 (/x 2 /m 2 ) in Figure |l(i|. The figure 
indicates the form factors are the same in the large /Lt 2 
limit. Since the leading-order contributions are identical, 
the Figure additionally shows the ratio of non-valence 
contributions to the form factor, namely 



5F I /8F J 



NI- 



(65) 



This ratio too tends to one as /1 2 becomes large, of course 
not as rapidly. Thus for finite [i 2 , we must add the cor- 
rection V — V in Eq. ( |45| ) which results in higher Fock 
contributions present in the non-instantaneous case. 

Lastly we remark that all of this is analytically clear: 
having eliminated different /1 2 dependence hidden in ip 
and </>, expressions for form factors in the instantaneous 
and non-instantaneous cases are identical to leading order 
in l//i 2 , i.e. 5J^,5Jg ,SJ^. ~ l//i 4 while the remaining 
non- valence contributions 5Jp,SJg match up with 8J^~ 
and 8J2 , respectively, to 1/fi 2 . 

Of course we have demonstrated the replacement of 
non-wave function vertices only to first order in pertur- 
bation theory. Schematically, we have required 



\V-V \ « \V(G-G)V\. 



(66) 
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FIG. 10: Comparison of form factors in instantaneous and 
non-instantaneous cases. The ratio Eq. d&l| ) of form factors 
(LO + NV) is plotted versus lqg 10 (^ 2 /m 2 ) at fixed £ = 0.707. 
Additionally the ratio Eq. ( |65| ) of non-valence contributions 
to the form factor (NV) is plotted. For fj, 2 > m 2 ,M 2 , the 
non-instantaneous contributions becomes approximately in- 
stantaneous. 



This is merely the condition that corrections to V from 
finite ji 2 are larger than second-order corrections from 
the reduction scheme in Eq. (phf). Quantitatively this 
condition translates to m/47r/i -C oe. If this condi- 
tion is met, the leading corrections are all of the form 
V — V and the instantaneous model becomes the ladder 
model (and non-wave function vertices disappear when 
calculating the form factor). When the condition is not 
met, but holds to [V(G — G)] 2 !/, say, we must work 
up to second order in the reduction scheme and use 
the second Born approximation. We do not pursue this 
lengthy endeavor here, since for the non-instantaneous 
case 7(x|M 2 ) oc 6[x(l—x)} necessarily excludes non-wave 
function contributions in time-ordered perturbation the- 
ory. In fact at any order in time-ordered perturbation 
theory we expect to find a complete expression for the 
form factor in terms of Fock component overlaps devoid 
of non-wave function vertices, crossed interactions, etc. 



V. SUMMARY 



these residues and postulate a form for the wave function 
are assuming not only a pole-free vertex, but a vertex in- 
dependent of light-front energy. These assumptions are 
very restrictive. 

This leads us formally to investigate instantaneous 
and non-instantaneous contributions to wave functions 
and form factors necessitating the reduction formalism 
Eq. J2q), We provide an intuitive interpretation for the 
light- front reduction, namely it is a procedure for approx- 
imating the poles of the Bethe-Salpeter vertex function. 
This procedure consists of covariant iterations followed 
by an instantaneous approximation cf Eq. (pq), where 
the auxiliary Green's function G enables the instanta- 
neous approximation, see Eq. (|2fj|). In order to calculate 
form factors in the reduction formalism, we construct the 
gauge invariant current Eq. (Bq). 



Using the ladder model (38) and an instantaneous ap- 
proximation (^4j) wc compare the calculation of wave 
functions and form factors in the light-front reduction 
scheme. Calculation of form factors is dissimilar for 
the two cases. In the ladder model, which is non- 
instantaneous, non-wave function vertices are excluded in 
time-ordered perturbation theory. For the instantaneous 
model, however, contribution from crossed interactions 
is required. Moreover, these instantaneous contributions 
derived are identical in form to non-wave function ver- 
tices used in constituent quark models. As a crucial check 
on our results, we take the limit fi 2 3> m 2 , M 2 for which 
the ladder model becomes approximately instantaneous. 
In this limit, calculation of the form factors for the two 
models is identical. 

The net result is an explicit proof, for (1 + 1) dimen- 
sions, that no n- wave function vertices are replaced by 
contributions form higher Fock states if the interactions 
between particles are not instantaneous. We call this a 
replacement theorem. The analysis relies on general fea- 
tures of the pole structure of the Bethe-Salpeter equation 
in light-front time-ordered perturbation theory. There- 
fore it is trivial to extend the theorem to (3 + 1) dimen- 
sions. The real question which remains is whether or 
not the interaction between light quarks can be approx- 
imated as instantaneous. We intend to learn how to use 
experimental data to answer this question. 



Above we investigate current matrix elements in the 
light-front Bethe-Salpeter formalism. First we present 
the issue of non-wave function vertices by taking up 
the common assumptions of light-front constituent quark 
models. By calculating the form factor in frames where 
£ + 1 0, Lorentz invariance mandates contribution from 
the vertex function's poles. Quark models which neglect 
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